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Abstract — We  investigate  optimal  estimation  for  both  the 
position  and  the  velocity  of  the  ground  moving  target  (GMT) 
by  employing  sensors  composed  of  unmanned  aerial  vehicles 
(UAVs).  The  problem  is  the  cooperative  sensing  by  the  UAVs, 
in  terms  of  their  location  geometries  to  achieve  optimal 
estimation  of  the  GMT.  Based  on  the  Cramer-Rao  bound, 
we  are  able  to  derive  the  minimum  achievable  error  variance 
in  estimation  of  the  position  and  the  velocity  of  the  GMT, 
and  obtain  the  optimal  geometries  of  the  UAV  sensors  via 
minimization  of  the  minimum  achievable  error  variance  for 
unbiased  estimation  commanded  by  the  Cramer-Rao  bound. 
Our  solution  is  complete  that  encompasses  various  situations 
for  the  GMT,  and  the  number  of  UAV  sensors. 

1.  Introduction 

Cooperative  sensing  with  UAVs  plays  an  important  role  in 
combat  ISR  (intelligence,  surveillance,  and  reconnaissance), 
such  as  GMTI  (ground  moving  target  identification)  and 
MSTE  (moving  surface  target  engagement).  The  objective 
is  the  early  identification  and  localization  of  the  threat  GMT 
of  interest,  estimate  and  track  the  GMTs  in  terms  of  their 
position  and  velocities.  The  ultimate  goal  is  the  destruction 
of  the  threat  GMT.  A  problem  arising  from  GMTI  and 
MSTE  is  the  sensor  geometries:  how  to  cooperate  the 
UAV  sensors  in  terms  of  their  locations  to  achieve  optimal 
estimation  of  the  position  and  velocity  of  the  GMT?  We 
assume  that  the  observed  data  by  UAV  sensors  are  radar 
signals,  comprised  of  noisy  azimuth,  range,  and  range 
rate  information.  For  simplicity  we  also  assume  that  the 
observed  data  are  corrupted  by  Gaussian  noise.  Based  on  the 
well-known  Cramer-Rao  bound,  we  provide  an  answer  to 
the  cooperative  and  optimal  sensing  problem.  We  will  show 
that  the  minimum  achievable  error  variance  for  estimation 
of  the  position  and  velocity  of  a  given  GMT  amounts  to 
maximization  of  certain  mathematical  quantity  that  involves 
the  angular  positions  of  the  UAV  sensors,  assuming  that 
all  UAV  sensors  have  the  same  distance  to  the  GMT.  The 
orthogonal,  or  near  orthogonal  sensor  geometries  for  UAVs 
are  preferred  in  most  of  the  applications.  We  will  provide  a 
complete  solution  to  the  optimal  sensor  geometry  problem 
with  an  arbitrary  number  of  UAV  sensors  that  encompasses 
various  situations  for  the  GMT. 

The  work  to  be  reported  in  this  paper  is  related  to  the 
general  problem  of  target  localization  [1],  [2],  [3],  [6],  [7], 
The  emphasis  of  this  paper  is  placed  on  the  geometries 
of  the  UAV  sensors  in  optimal  estimation  of  not  only 


position,  i.e.,  localization,  but  also  velocity  of  the  GMT 
due  to  tracking,  which  is  relatively  new.  We  consider  the 
case  of  two-dimensional  plane  for  the  coverage  field,  which 
is  adequate,  if  the  elevation  of  the  field  is  known  a  priori, 
that  is  true  in  most  of  the  applications.  The  problem  of 
cooperative  sensing  is  also  related  to  collaborative  signal 
processing  [4]  in  the  sense  that  sensors  need  work  together 
to  achieve  optimal  estimation.  The  notations  in  this  paper 
are  standard.  The  inner-product  and  the  Euclidean  norm 
of  column  vectors  are  denoted  by  <  >.  and  ||  •  || , 

respectively,  and  the  transpose  by  -T.  Other  notations  will 
be  made  clear  as  we  proceed. 


2.  Cooperative  Estimation 


Let  the  number  of  UAV  sensors  be  m  >  1,  with  their 
respective  positions  and  velocities  known.  For  the  fcth  UAV, 
its  position  and  velocity  at  time  t  are  denoted  by  pk(t),  and 
vk(t),  respectively.  Let  pr(t),  and  vrit)  be  the  position  and 
velocity  of  the  GMT  respectively.  Then 

Prif)  =  Pk  {t)+6pk{t),  vT(t)  =  =  vk(t)+6vk(t ), 

at 

with  dpk(t)  and  5vk{t)  being  the  respective  differences  of 
the  position  and  velocity,  respectively,  between  the  GMT 
and  the  fcth  UAV  Let  the  rc-axis  be  associated  with  east- 
west,  and  y- axis  associated  with  north-south.  Then  we  have 
the  following  representations: 


Pk(t )  = 

Spk(t)  = 


Xk(t) 

.  Vk(t)  \  ’ 
6xk(t) 

Syk(t)  J  ’ 


Vk(t)  = 

Svk(t)  = 


Xk(t) 
Vk{t)  \  1 
Sik(t) 
Syk{t) 


Let  9k(t)  =  L 8pk(t )  be  the  argument  of  8pk{t),  that  is  the 
bearing  parameter  between  the  kth  UAV  and  the  GMT  (cf. 
Fig.  1). 


Fig.  1  The  GMT  and  the  observable  cone  by  UAV 
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The  radial  velocity  of  the  GMT  observed  by  the  fcth  UAV 
induces  the  Doppler  shift,  denoted  by  u>k(t),  which  is  given 
by 

Wfc(t)  =  cos  (Ok(t))xT(t)  +  sm(0k{t))yT{t).  (1) 

For  simplicity  we  will  skip  the  time  variable  t  in  sequel. 
We  comment  that  u>k  is  the  orthogonal  projection  of  vr 
onto  5pk,  i.e.,  u>k  =  |Ft||  cos  (9k  —  Ivt)-  The  greater  the 
velocity  vt,  the  easier  for  radars  to  capture  its  movement 
in  measuring  its  azimuth,  range,  and  range  rate.  Hence  the 
observation  of  the  range  rate  by  the  fcth  UAV  sensor  has 
the  form  of 


Zk(oj)  =  /(wfe)wfc  +  T]k(u>),  (2) 


where  rjk(u>)  is  the  Gaussian  noise  with  variance  <r^  , 
and  /(•)  is  a  nonlinear  function  characterizing  the  GMT 
indication.  As  a  rule  of  thumb,  the  GMT  is  undetectable  if 
ujk  is  smaller  than  10  kilometers  per  hour,  which  is  taken 
as  the  critical  velocity,  denoted  by  u>tc- 

Suppose  that  u)k  is  large  enough  such  that  the  GMT  is 
observable,  or  f(cuk)  =  1.  Then 


zk{x) 

.  zk{y)  _ 


XT 

Vt 


Vk(x) 

.  Vk(y)  _ 


(3) 


are  observations  of  px  by  the  fcth  UAV  sensor  with  r)k{x) 
and  t) k{y)  white  Gaussian  noise  of  zero  mean.  If  Ok  =  0, 
then  rjk  ( x )  and  rjk  (y)  are  the  corruption  noises  in  measuring 
the  range  and  azimuth  respectively.  Thus  ijk(x)  and  ijk(y) 
are  uncorrelated  with  variances  cr^fe  =  apk  (range  variance) 
and  o-yk  =  a^k  (azimuth  variance),  respectively.  Or  the  joint 
probability  density  (JPD)  is 


P[Vk(x),Vk(y)] 


- - exp  I  —  - 

27T(Jr,  (Jn,  \  2 


vl(x) 


That  is,  the  constant  probability  density  contour  has  an 
ellipsoidal  shape.  But  if  Ok  ^  0,  then  rjk (x)  and  rjk{y)  are 
correlated,  and  the  covariance  matrix  is  given  by 


^ek  =  U(0k) 


rk 

0 


U(0k)T, 


where  (/(•)  is  a  rotation  matrix,  given  by 


U(0) 


COS  (0)  —  8111(0) 

sin(0)  cos(0) 


It  follows  that  the  the  IPD  has  the  form  of 


exp 


P[Vk(x),Vk{y)]  = 


Vk{x) 

Vk(y) 


-,T 


^-1 


Ok{x) 

Vk(y) 


27rx/det(S^J 


(4) 


It  should  be  clear  that  the  unitary  matrix  U(0)  performs 
rotation  of  angle  0.  For  simplicity  we  assume  that  each 
UAV  sensor  has  the  same  noise  variances  for  measuring  the 
azimuth,  range,  and  range  rate,  respectively,  if  the  distances 


to  the  GMT  are  the  same,  which  are  true,  provided  that 
identical  sensors  are  used  in  each  UAV. 

The  position  measurement  (3)  assumes  that  ojk  >  u>c.  In 
the  general  case,  the  observation  of  the  position  of  the  GMT 
in  (3)  needs  be  replaced  by 


Zk{x) 

.  zk{y)  _ 

£1  E-h 

1 _ 1 

II 

/(wfc)  + 

Vk  (x) 

_  Vk{y)  . 

where  the  noise  has  the  JPD  as  in  (4). 

Denote 

r 

.  iT 

X  =  [  XT  yr  XT  Vt  \ 

'  0  0 

cos(0fc) 

sin (0k)  ' 

C’k  =  f{uk) 

1  0 

0 

0 

0  1 

0 

0 

’  ' 

’  zk(uj)  ' 

%  = 

Vk(x) 

5  —k 

Zk(x) 

.  Vk(y)  . 

.  zk(y)  _ 

Combining  (2)  and  (5)  yields 


(5) 

(6) 

(7) 

(8) 


Zk  =  ckX  +  Vk,  ^k  =  ^iUkVl}  =  dias{^lk,^ek), 

(9) 

for  1  <  fc  <  m.  Because  /(•)  and  0k  are  nonlinear  functions 
of  the  parameter  vector  \  t0  be  estimated,  the  observed 
data  vector  zk  is  also  a  nonlinear  function  of  the  parameter 
vector  x  in  general.  However  it  is  a  special  nonlinear 
function  of  the  parameter  vector.  Indeed  if  0k  is  known, 
and  l Ok  >  <W,  then  the  observed  data  vector  zk  becomes  a 
linear  function  of  the  parameter  vector  %.  It  is  this  quasi¬ 
linearity,  which  will  be  exploited  in  this  paper  to  derive  the 
Cramer-Rao  bound  for  the  underlying  estimation  problem. 

It  is  noted  that  Ck  is  a  function  of  Ok-  Hence  if  all  UAVs 
have  the  equal  distance  to  the  GMT,  then  the  geometries  of 
the  m  UAVs  are  uniquely  determined  by  {Ok}™^-  Recall 
that  Ok  is  the  bearing  parameter  from  the  fcth  UAV  to 
the  GMT.  Our  cooperative  sensing  problem  is  to  find  the 
optimal  geometries  of  m  UAVs  so  that  the  minimal  error 
variance  of  the  estimation  error  is  minimized.  Denote 


= 

'  Z.1 

,  c  = 

1 

...  cy 

_ 1 

,  A/-  = 

'  Hi  ' 

.  —m  . 

1 

P 

1 _ 

-  ILm  - 

It  follows  that 

tlr  =  CX+  M.  (10) 


In  light  of  (9),  the  Gaussian  assumption  implies  that  ip  has 
a  JPD  function 


eXP  \  te*  -  °kX)T  1  {zk 


fc= 1 


(27r)2m  det(Efc)  2 


The  error  covariance  X  corresponding  to  %,  any  unbiased 
estimate  of  X,  based  on  observations  ip  =  CX  +  A f, 
is  lower  bounded  by  the  well  known  Cramer-Rao  bound 
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matrix  Aqr;  Its  inverse  XfA\,  is  the  Fisher  information 
matrix.  That  is,  A  >  Acr  for  any  unbiased  estimation 
algorithm.  Furthermore  the  Gaussian  distribution  implies 
that  the  Cramer-Rao  bound  is  achievable  asymptotically. 
Our  motivation  is  to  minimize  Acr  by  choosing  { 6 &} 
correctly,  which  will  yield  the  optimal  geometries  of  the 
UAVs  under  the  assumption  of  fixed  equal  distance  between 
the  GMT  and  each  UAV.  However  minimization  of  Acr  is 
a  daunting  task.  Instead  we  opt  to  minimize  Tr{  Acr  }.  that 
is  more  feasible.  We  notice  that  for  any  unbiased  estimate 
X,  there  holds 

E{||X-X||2}  =  Ti-{E[(X-x)(x-x)T]} 

=  Tr{A}  >  Tr{ACR}. 

Consequently  minimization  of  TrjAcR}  makes  more  sense 
in  searching  for  the  optimal  geometries  of  the  UAVs,  which 
will  be  studied  in  later  sections. 


3.  The  Fisher  Information  Matrix 

In  this  section  we  compute  the  Fisher  information  matrix, 
associated  with  Cramer-Rao  bound  based  on  observation 
ip  =  Cx  +  A/",  by  assuming  that  the  distances  from  UAVs  to 
the  GMT  are  all  the  same,  and  u>k  is  large  enough  such  that 
each  UAV  is  inside  the  observable  cone  for  which  f(oJk)  ^ 
0,  and  thus  Cf  ^  0  for  1  <  k  <  m.  Let  n  (=  4)  be  the 
dimension  of  x-  According  to  the  Slepian-Bangs  formula 
[7]  (page  289),  the  Fisher  information  matrix  is  given  by 


X 


-l 

CR 


-E 


9%dxT  J  L2  l  dXj  J 


+ 


dCx\Ts-i  fd£x\ 

dXi  )  V  dXi  J 


where  E  =  diag(Ei,  E2,  •  •  • ,  Em),  and  Xfc  is  the  fcth  element 
of  x-  Because  of  the  diagonal  structure,  we  can  compute  each 
diagonal  block  separately. 

The  equal  distance  assumption  for  all  the  UAV  to  the  GMT 
implies  that  ||5p*;||  =  ||<5p||,  and 


in  light  of  tan(f?)  =  ( yr  —  y)/(xr  —  x).  Denote 


Tg  = 


0 


U{9)rd^-U{9) 


0 


0 

.-1 


By  the  property  of  trace,  and  the  expression  of  Eg,  we  have,  in 
the  case  of  multiple  UAV, 


Tr 


Z>{rU 

11=1 


ddfi  d6fi 

dXi  dxj  ' 


We  choose  to  skip  the  complex  expression  of  Tr  {T2},  and  of  the 
derivatives  of  ip  with  respect  to  Xi- 

We  notice  that  by  lot  =  u>t{9)  =  X3  cos (9)  +  X4  sin (9), 

[X4  cos(6»)  -  X3  sin(0)]2  +  [x3  cos (9)  +  X4  sin(0)]*  =  xi  +  xl- 


It  is  appropriate  to  define  u>t±  ( 9 )  =  \4  cos (9)  —  X3  sin(0), 
which  is  the  orthogonal  projection  of  vt  to  the  vector  that  is 
perpendicular  to  Sp,  while  lot  is  the  orthogonal  projection  of  Vt 
to  the  vector  Sp.  Hence  lot(6)2  +lot±(9)2  =  llnrll2.  Let 


9 


(fc)  = 


-  sin(6L) 
cos (9k) 


9 W  = 


cos (9k) 
sin  (9  k) 


Then  we  obtain  the  expression  for  the  Fisher  information  matrix 
associated  with  the  Cramer-Rao  bound: 


AT 


2||  Sp\[- 


■E 

k=i 


Tr{Lifc} 


Q(k) 


0 


+ 


E 

k=l 


0 


+E 


lot ( 9k)9 . 


(fc) 


LOT 


{9k)9<£> 

»p||  9(£> 


Remark  3.1:  It  is  worth  to  pointing  out  that  lot±  <  ||'Ut||, 
which  is  the  driving  speed  of  the  GMT.  Hence  lot±  «  ||<5p]| 
holds.  For  this  reason,  the  Fisher  information  matrix  correspond¬ 
ing  to  the  Cramer-Rao  bound  has  an  approximate  block  diagonal 
form  as  A'cr  ~  diag  ^A'^,  A^^  where  A^  and  A^  are 
the  Cramer-Rao  bound  matrices  for  estimation  of  velocity  and 
position,  respectively.  That  is,  the  Cramer-Rao  bound  matrices  for 
velocity  and  position  are  decoupled,  and  given  respectively  by 


LOk  =  LOT, 


(7  ru  —  O' r 


where  1  <  k  <  m.  For  the  case  of  single  UAV  (m  =  1),  E  = 
diag (cr 2 ,  Eg),  and 

Eg  =  U(9)dmg(o2r,o2a)U(9)T,  (11) 

Eg"1  =  U(9)dmS(af2,of2)U(9)T.  (12) 

Direct  computation  gives 

dEg  /  2  2\ 

and  =  0  V  i  in  which  we  have  assumed  that  crpf  is  independent 
of  x-  Since  9  is  a  function  of  (xi>X2)  =  ( xt,TJt ),  but  not  a 
function  of  (X3,X4)  =  ( xt,Vt ), 


VVR  = 

Er.A.-,r‘ 


A, 


(p) 

CR 


The  approximate  expression  for  A'q1^  is  due  to  ||5p||  >>  Tr  {Tgfc  } 
for  1  <  k  <  m.  I 

In  the  following  we  consider  the  Cramer-Rao  bounds  for 
estimation  of  position  and  velocity  separately,  in  light  of  Remark 
3.1.  For  velocity  only,  the  Fisher  information  matrix  has  an  inverse 


with 


sin(2f?) 

—  cos(2 9) 

as 

—  cos(2 9) 

sin(20) 

5 

v(v)  _  2 

ACR  — 

||c||2  <C,S> 

_  <  c,s  >  ||s||2 

where  s  =  [  si  •  •  •  sm  \ T  and  c  =  [  ci  •  •  •  Cm 
Sk  =  sin(#fc)  and  ck  =  cos(6L).  It  can  be  verified  that 


( Vt  -  y) 


sin  (9) 


d9 

dx  1  (xt  -  x)2  +  (i/t  -  y)2 
d9  _  (xt  —  x)  _  cos(6>) 

dx 2  (xt  -  x)2  +  (yT  -  y)2 


d9 

dx  3 

d9 

dx  4 


=  0, 
=  0, 


Tr 


(x&) 


IUII2-|<c,s> 


(13) 


Clearly  the  above  is  minimized,  if  and  only  if  its  denominator 
is  maximized  over  which  determine  the  geometries  of 
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the  m  UAV  sensors.  For  position  estimation  only,  the  minimal 
achievable  error  variance  based  on  the  Cramer-Rao  bound  is  given 
approximately  by 


where  with  So2  =  (ay,  2  —  oa  2)2,  den  is  given  by 


4.  Optimal  UAV  Sensing  Geometries 

As  shown  in  the  previous  section  that  the  minimization  of  the 
minimum  achievable  error  variance  among  all  unbiased  estimation 
algorithms  is  equivalent  to  maximization  of  (||c||2||s||2  — |  <  c,s> 
|2)  over  the  UAV  sensor  geometries  {9k}-  This  problem  will  be 
solved  in  this  section.  For  convenience,  we  denote  A r-y  =  9k  —9i, 
and  A  as  the  angle  of  the  observable  cone.  There  holds 


den  =  m2 ay.  2aa  2  +  8a2  (||c|j2||s||2  —  |  <  c,  s  >  |2)  . 

We  come  across  with  (||c||2||s||2  —  <  c,s  >  |2)  once  more. 

Flence  for  m  >  2, 

Tr{XCn}  =  TV  {xg }  +  Tr  {*<£>}  (15) 

mo2  m(oyT2  +  a“2) 

||c||2||s||2  —  |  <  c,  s  >  |2  den 

The  approximation  is  due  to  the  large  value  of  the  distance  from 
the  UAV  sensors  to  the  GMT,  which  holds  true.  The  next  result 
follows. 

Proposition  3.2:  Suppose  that  the  m  UAVs  have  equal  distance 
||<5p||  to  the  GMT  with  ||<5p||  sufficiently  large  compared  with 
azimuth  and  range  variances,  and  the  velocity  of  the  GMT.  Let 
j,  which  specify  the  geometries  of  the  m  UAV  sensors. 
Then  minimization  of  the  achievable  minimum  error  variance 
Tr  {Xch}  over  all  possible  unbiased  estimation  algorithms  is 
equivalent  to  maximization  of  (||c||2||s||2  —  |  <c,s>  |2)  over 
the  sensor  geometries  {9k  Yk=i-  provided  that  all  m  UAV  sensors 
are  within  the  observable  cone. 

The  above  proposition  states  that  optimal  cooperative  sensing 
is  equivalent  to  the  mathematical  problem  of  maximization  of 
( II c|| 2 1 1 s || 2  —  |  <c,s>  | 2 ),  which  will  be  investigated  in  the  next 
section.  Before  ending  this  section,  we  would  like  to  study  this 
maximization  issue  for  small  values  of  m  as  follows.  In  the  case 
of  m  —  2, 


llcflkll2  -  I  <  c,s  >  |2  =  sin2(6>i  -  02), 


l[c||2||s||2-|<c,s>|2=  ^2  sin2(A fe,i).  (16) 

k>i= 1 


Based  on  the  analysis  in  the  previous  section  for  m  =  2,  and 
m  =  3,  it  is  expected  that  the  optimal  UAV  sensing  geometries 
are  easier  to  determine,  if  the  observable  cone  has  an  angle  A  no 
greater  than  90°.  This  is  indeed  true,  as  stated  in  the  next  theorem. 

Theorem  4.1:  Denote  gm(9i,  d2,  •  •  • ,  9m)  =  ||c||2||s||2  —  |  < 
c,s>  j2.  For  A  <  7t/2,  we  have 


s  max  gm(9i,- ■  ■  ,9m)  = 


12^— t  sin2 (A),  odd  m 
yj-sin2(A),  even  m 

(17) 


Proof:  Consider  first  the  odd  m.  If  we  take  A k,i  =  A  for 
1  <  i  <  p,  and  p  +  1  <  k  <  m  where  p  =  (m  —  l)/2  or 
p  =  ( m  +  l)/2,  then 


Qm{9 1 , 02,  '  '  '  ,  9m'} 


sin2  (A). 


It  is  claimed  that  the  above  is  the  maximum  of  gm{-,  We 

use  induction.  The  claim  is  true  for  m  —  3,  as  shown  in  the 
previous  section.  For  m  >  3  odd. 


Qm  {9\ ,  *  *  •  ,  9m)  —  9m  —  2(^2,  '  "  ,  9m  —  l)  T  h{9\,  •  •  •  ,  9m} 
where  h{6 1, 92,  ■  ■  • ,  9m)  has  the  expression 

m  —  1 

h{-  ■  ■}  =  sin2(Am>1)  +  y^[sin2(Amy)  +  sin2(Ai,i)].  (18) 

i= 2 

For  1  <  i  <  m,  A m,i  +  Ai}1  =  Amjl  and  A m,i  —  A,,!  = 
Am>1  -  2Ai,r.  Thus  with  d<(Ai,i)  =  sin2(Amy)  +  sin2(A;,i), 


leading  to  the  orthogonal  geometry  (9\  —92)  =  90°  easily  for  the 
optimal  solution.  Thus 

Tr  {Xcr)  >  2o-2  H - — ^2  • 

(Tr  4 ~  tT a 

The  right  hand  side  is  attained,  if  the  observable  cone  has  an  angle 
no  smaller  than  90°. 

In  the  case  of  m  =  3,  it  can  be  shown  that 


di(Ai,i)  =  1  -  cos(Am,i)  cos(Am:i  -  2Aiji). 

Without  loss  of  generality  Amii  >  0  is  assumed,  which  is 
equivalent  to  9m  >  9\.  Since  Am>i  <  A  <  7r/2,  di(Ayi)  is 
maximized  by  taking  Ayr  =  0  and  Am,i  =  A,  leading  to 

maxdi(Ai  1)  =  sin2(Am  1)  <  sin2(A) 

Ayr 

=>  max  h(6i,- ■  ■  ,9m}  =  (m  —  1)  sin2(A). 


max{||c||2||s||2  -  |  <  c,s  >  |2} 

Ok 

{2.25,  A  >  120°, 

2 sin2 (A),  A  <90°, 

2  +  2  cos(A)  cos2(A/2),  90°  <  A  <  120° 

where  A  is  the  angle  of  the  observable  cone.  In  the  general  case 
of  m  >  2,  we  have 

m  k  —  1 

l|c||2||s||2  -  |  <C,S>  I2  =  ^^sin2(6ifc  -9i). 

k=2  i=  1 

For  optimal  cooperative  sensing,  we  need  maximize  the  above 
summation,  which  will  be  studied  in  next  section. 


By  induction  assume  that  gm-2{92,  ^3,  •  •  • ,  9m—  1)  is  maximized 
by  taking  Ayi,  =  A  for  2  <  i  <  p,  and  /r+l<fc<m—  1 
where  p  =  ( m  —  l)/2  or  p  =  ( m  +  l)/2  with  maximum  [(m  — 
2)2  —  1]  sin2(A)/4.  With  9\  and  9m  added  in,  we  obtain 


maxgm(0r,02,  •  •  •  ,6m)  <  maxgm_2(02,  •  •  • ,  9m-i} 

+  max  h(9i,  92,  ■  ■  ■  ,9m} 


m  —  1  .2 


sin"  (A), 


which  is  indeed  achieved.  For  the  case  of  even  m,  we  assign 
{QkYkiX  on  one  edge  of  the  cone,  and  +1  on  the  other 

edge  of  the  cone,  leading  to 


gm{9i,92,-  ■  ■  ,9m}  =  sin2  (A)  =  ^-sin2(A). 
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The  same  induction  procedure  as  earlier  can  be  used  to  prove  that 
the  above  is  the  maximum,  thereby  proving  (17)  for  the  case  of 
even  m,  which  is  omitted.  B 


Theorem  4.1  favors  clearly  the  orthogonal  assignment  for 
optimal  geometries  of  the  UAV  sensors,  if  the  observable  cone 
has  an  angle  no  smaller  than  90°.  Otherwise  it  favors  to  assign 
roughly  half  of  the  UAV  sensors  on  one  edge  of  the  cone,  and  the 
other  half  or  so  on  the  other  edge  of  the  cone.  While  the  result 
in  Theorem  4.1  is  intuitively  appealing,  the  induction  procedure 
in  the  proof  is  not  extendible  to  the  case  of  A  >  90°.  A  different 
solution  approach  is  necessary. 

In  light  of  Theorem  4.1,  m2/ 4  seems  to  be  an  upper  bound  for 
gm{9 1,  •  •  •  ,8m),  if  A  >  90°,  which  is  true  for  the  case  of  even 
m.  The  next  result  is  more  general. 

Lemma  4.2:  Suppose  that  6k  E  [0,  n)  are  arbitrary  for  fc  = 
1,  2,  •  •  • ,  n.  Then  there  holds 

m  2 

Am  (01, 02,  ■  ■  ■  ,  dm)  =  sin2(0fc  —  9i)  < 

k>i=  1 

In  addition  the  above  upper  bound  inequality  is  tight. 

Proof:  Let  <j>k  =  2 6k  for  all  k.  Then 

m  m 

2gm{9i,  02,  •  •  • , 6m)  =  ^^sin2^-^)  (19) 

k= 1  i= 1 

2  mm 

=  it  _  ^EEcos^fc 

k= 1  i= 1 

It  is  claimed  that  the  following  lower  bound 

n  n 

EE  cos((pk  -  (pi)  >  0  (20) 

fc=l  i  =  l 

holds,  and  is  tight.  Indeed  we  have 


0  < 


n  k  —  1 


n  + 


k= 2  £=1 


n  k  —  1 

n  +  2  EE  cos((pk  -  (pi) 

k= 2  i=  1 


n  n 

=  EE  COS (<f)k  -  (f)i) 

k= 1  i= 1 


which  concludes  the  proof  of  (20).  Furthermore  by  taking  (pk  = 
2kn/n  for  all  k ,  we  have 


EE  cos(0jt  -  (j)i)  =  0, 

k= 1  i= 1 

that  shows  that  the  lower  bound  in  (20)  is  tight.  It  follows  from 
(19)  that  gm(0i,02,  •  •  •  ,6m)  <  m2/ 4,  which  is  tight,  by  taking 
6k  =  kn/n,  in  light  of  (20)  again.  0 

The  upper  bound  established  in  Lemma  4.2  assumes  that  {6k  } 
are  arbitrary,  considering  that  sin2(-)  is  periodical  with  period  i r. 
In  our  case,  |Afc,j|  =  |0fc  — 0;|  are  restricted  to  be  within  A,  which 
can  be  strictly  and  considerably  smaller  than  n.  Hence  m2/ 4  may 
not  be  achievable  as  shown  in  Theorem  4.1  in  the  case  of  odd  m, 
even  if  A  =  n/2,  that  is  rather  different.  However  in  the  case  of 


even  m,  we  obtain  the  following  rather  easily  in  light  of  Lemma 
4.2. 

Corollary  4.3:  For  A  >  n/2  and  even  m, 

flm(01,02,---,0m)  =  ||c||2||s||2  ~  |  <  C,  S  >  |2  <  m?  /  4, 
which  is  attained. 

The  above  result  follows  easily  from  the  fact  that  m2/ 4  is  an 
upper  bound,  that  is  achieved  for  the  case  A  =  n/2  and  even  m. 

Remark  4.4:  There  are  more  than  one  assignment  of  {0fc}  in 
achieving  the  upper  bound  m2/ 4,  if  A  >  n/2  and  m  is  even. 
Indeed  by  grouping  {Sfcjj’Aj  into  m/2  pairs,  with  each  pair- 
satisfying  A k,i  =  n/2,  the  upper  bound  m? / 4  can  be  achieved, 
regardless  of  the  location  of  each  pair  so  long  as  they  are  within 
the  observable  cone,  which  is  illustrated  in  Fig.  2  with  m  =  4.  This 
fact  has  its  root  in  the  necessary  condition  of  optimality.  Setting 
its  partial  derivatives  with  respective  to  {6k}  to  zero  yields 


<9ffm 

ddk 
dgm 
36 1 


3gm— 2 

36k 


2  cos(Am,i)  sin(Am,fc  -  Afc,i) 


sin(2Am,i) 


m  —  1 

1  ^  ^  COS (Am,k  ^k, l) 

k= 2 


dgm 

36m 


m  —  1 

cos(Am,i)  sin(Am,fc  -  Afc,i) 

k= 2 


0, 

0, 

0. 


Hence  with  optimally  assigned,  then  the  enforcement  of 

Am>i  =  dm  —0i  =  90°  will  satisfy  the  above  necessary  condition, 
regardless  of  the  values  of  d\  and  dm,  that  happen  to  be  globally 
optimal  for  even  m  when  A  >  n/2.  This  is  advantageous  in 
applications,  because  the  geometries  of  the  first  (m  —  2)  UAVs 
need  not  be  changed,  if  an  additional  pair  of  UAVs  is  added  in.  I 


Vj 


The  difficult  case  is  clearly  odd  m,  for  which 
gm(9 1, 02,  •  •  • ,  dm)  has  the  maximum  (m2  —  l)/4,  if  A  =  n/2 
as  shown  in  Theorem  4.1.  On  the  other  hand,  when  A  is  close 
to  7r,  gm (0i, 02,  ■■■  ,9m)  has  m2/ 4  as  the  maximum,  by  Lemma 
4.2.  The  following  result  is  constructive. 

Theorem  4.5:  Let  gm(di,  02,  •  •  •  ,  9m)  be  as  in  Theorem  4.1 
with  m  odd.  If  A  >  2n /2>,  then 

max  pm (0i ,  02,  •  •  • ,  dm)  =  m2/4. 

^1  >  ■  ■  ‘  )0m 

Proof:  We  need  only  show  that  m2/ 4  is  achievable,  because  it 
is  an  upper  bound  by  Lemma  4.2.  For  m  =  3,  the  result  from  the 
previous  section  can  be  used  to  show  that 


m= 3 

max  sin2(0/c  —  Of)  =  3sin2(7i/3) 

01,02,0+3  ' 

k>i= 1 


3^ 

4  ’ 


which  holds  by  taking  A3^  =  A24  =  n/3  and  A3.1  =  2n/3. 
For  m  >  3  odd,  we  set  {0;}^.!  as  illustrated  in  Fig.  3  with 
n  =  (m  —  l)/2.  It  is  claimed  that  gm  (0i ,  •  •  • ,  0m)  =  m2/ 4. 
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Using  the  induction  consider  first  m  =  5.  We  have 

m= 5 

sin2(0fc  —  6i)  =  6.25. 

k>i= 1 

Thus  the  result  is  true  for  the  case  m  =  5.  For  m  >  5  that  is 
odd.  the  decomposition  in  the  proof  of  Theorem  4.1  can  be  used 
to  arrive  at 

gm(Ol,0  =  (m~2)2  +  h(0l,92,  •  ■  •  ,  6m), 

where  gm-2(92,  ■  ■  ■  ,  8m-i)  =  (m  —  2)2/4  is  assumed  by  the 
induction  argument,  and 

h(9 1,02,  •  •  • ,  0m)  =  (m  -  1), 


by  cos(2(^0)  =  — 1  /(m  —  1).  11 

It  is  noted  that  2<^0  =  2-7r/3  =  120°  for  m  =  3,  that  coincide 
with  the  result  in  the  previous  section,  and  that  in  Theorem  4.5. 
For  m  =  5,  we  have  2ip0  =  104.5°.  Hence  as  m  increases,  2ipa 
approaches  to  n/2  =  90°,  consistent  with  the  result  for  even  m. 
The  next  result  deals  with  the  case  of  odd  m  with  90°  <  A  <  2p0, 
which  is  the  last  scenario  uncovered. 

Corollary  4.7:  Suppose  that  m  is  odd,  and  ip0  €  \ir  / 4,  n/2) 
satisfies  1  +  (m  —  1)  cos(2ip0)  =  0.  If  90°  <  A  <  2 pa  <  120°, 
then  with  n  =  ( m  —  1)  there  holds 

max  ■  ■  ,9m)  =  j  [4sin2(A/2)  4-  nsin2(A)l  . 

The  proof  of  this  corollary  is  skipped  due  to  page  limit. 


by  A m,i  +  A;,!  =  Am,i  =  7r/2.  It  follows  that 


gm  {9l ,  92  ,  '  '  '  j  9lla)  — 

that  concludes  the  proof. 


(m  -  2)2 


0|i+2  A 

0114-S 


+  (m  -  1)  =  — 


^  0n 

01,...,n-l 


Fig.  3  Optimum  sensing  geometry  for  m  =  2n  +  1 

Fig.  3  shows  that  there  are  (n  —  1)  =  (m  —  3)/2  orthogonal 
pairs.  In  fact  these  ( n  —  1)  orthogonal  pairs  need  not  be  in  the 
same  location  as  illustrated  in  Fig.  3  for  which  the  optimality  still 
holds,  that  follows  from  the  same  argument  as  in  Remark  4.4.  For 
the  case  7t/2  <  A  <  2n/Z  and  odd  m,  the  situation  is  more 
delicate.  Nevertheless  we  have  the  following  result. 

Theorem  4.6:  Suppose  that  m  is  odd  and  ipa  £  [tt /4,  n/2) 
satisfies  1  +  (m  —  l)cos(2<^0)  =  0.  Then  for  A  >  A0  =  2i p0, 
there  holds 

max  gm(0i,  92,  ■  ■  ■ ,  9m)  =  m2/4 

Proof:  We  choose  {9k}™=i  such  that 
02e  —  92i-i  =  2  <p,  d2e  —  dm  =  9m  —  92t-\  =  v  (21) 
for  i  =  1, 2,  •  •  • ,  (m  —  l)/2.  Then  it  can  be  verified  that 

gm(9 1,  •  •  ■ ,  9m)  =  (y,n  [4 sin2 ((^5)  +  (m  -  1)  sin2(2^)]  , 
which  is  denoted  by  w(ip).  Setting 

=  (m  —  1)  fsin(2(p)  +  (m  —  1)  sin(2y>)  cos(2y>)]  =  0 

dip 

yields  sin(2^)  =  0  or  l  +  (m—  1)  cos(2y>)  =  0.  We  are  interested 
in  the  case  0  ^  2ip  <  7r,  by  the  fact  that  2ip  <  A  <  n.  Hence 
ip  =  ip0  is  the  maximal,  yielding 

w(<A>)  _  -„2/  x  .  (m-  1\  2M,  s 


in2(v^c)  +  I  1  j  sin2 (2^0 ) 

L  -costMy  +  (m^l)  (l_cosJ(2„o)) 


1  —  cos(2y>0)  i  ( m  —  1 

2  +  V  4 

,  .  m2 

w(<Po)  = 


5.  Conclusion 

We  studied  the  cooperative  sensing  problem  for  estimation  of 
the  position  and  velocity  of  the  GMT  by  employing  multiple  UAV 
sensors.  The  measured  data  are  radar  signals  of  azimuth,  range,  and 
range  rate  corrupted  by  Gaussian  noise.  Assuming  equal  distance 
from  UAV  sensors  to  the  GMT,  the  geometries  of  the  UAVs  are 
uniquely  determined  by  their  angular  positions.  We  have  employed 
the  Cramer-Rao  bound  to  obtain  the  minimum  achievable  error 
variance  based  on  the  assumption  that  the  measurement  noise 
is  Gaussian  to  which  the  Slepian-Bangs  formula  applies.  The 
minimum  achievable  error  variance  is  a  function  of  the  UAV 
sensors  geometries,  that  favors  orthogonal,  or  near  orthogonal 
assignment  for  the  UAV  sensors.  We  comment  that  except  in  the 
case  of  odd  m  and  2ipa  <  A  <  120°,  the  optimal  UAV  geometries 
admit  the  property  that  no  change  needs  be  made  for  the  existing 
UAVs,  if  an  additional  pair  of  UAV  sensors  is  added  into  the 
existing  ones,  that  is  why  induction  argument  is  used  in  the  proof, 
as  commented  in  Remark  4.4.  This  is  advantageous  in  applications. 
However  in  the  case  of  odd  m  with  90°  <  A  <  120°,  this 
advantage  is  lost.  The  geometries  of  all  the  UAV  sensors  need 
be  reconfigured,  when  new  ones  are  added  in.  We  would  like  to 
point  out  that  the  optimal  cooperative  sensing  problem  does  not 
take  the  trajectories  of  the  UAV  sensors  into  account,  which  is 
currently  under  study. 
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